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The physics of manganites is often described within an effective two-band tight-binding (TB) 
model for the Mn e g electrons, which apart from the kinetic energy includes also a local "Hund's 
rule" coupling to the ti g core spin and a local coupling to the Jahn- Teller (JT) distortion of the 
oxygen octahedra. We test the validity of this model by comparing the energy dispersion calculated 
for the TB model with the full Kohn-Sham band-structure calculated within the local spin-density 
approximation (LSDA) to density functional theory. We analyze the effect of magnetic order, JT 
distortions, and "GdFe03-type" tilt-rotations of the oxygen octahedra. We show that the hopping 
amplitudes are independent of magnetic order and JT distortions, and that both effects can be 
described with a consistent set of model parameters if hopping between both nearest and next- 
nearest neighbors is taken into account. We determine a full set of model parameters from the density 
functional theory calculations, and we show that both JT distortions and Hund's rule coupling are 
required to obtain an insulating ground state within LSDA. Furthermore, our calculations show 
that the "GdFeOa-type" rotations of the oxygen octahedra lead to a substantial reduction of the 
hopping amplitudes but to no significant deviation from the simple TB model. 



I. INTRODUCTION 

Manganite systems, Ri- x A x MnOs, where R is a triva- 
lent rare earth cation (e.g. La 3+ , Pr 3+ , Nd 3+ , . . . ) and 
A is a divalent alkaline earth cation (e.g. Sr 2+ , Ca 2+ , 
. . . ) , have attracte d the attention of scientists already 
for decades? 1 * 2 * 3 * 4 * 5 ^ These compounds exhibit a very rich 
phase diagram as a function of both temperature and 
composition, with various types of eventually coexisting 
charge, orbital, and magnetic order, and they are there- 
fore important prototype materials to test our current 
understanding of correlated electron systems. In addi- 
tion, the observation of "colossal magneto-resistance" P 
a magnetic-field induced change in electric resistivity by 
several orders of magnitude, has spawned further interest 
both in the fundamental physics behind this effect as well 
as in the question of whether this effect can be utilized 
for technological applications. 

LaMnOa , the parent material for many manganite sys- 
tems, exhibits an orthorhombically distorted perovskite 
structure with Pnma space group (see Fig. [tk)P The 
observed deviation from the ideal cubic perovskite struc- 
ture (shown in Fig. [TJj) involves both Jahn- Teller (JT) 
distortions of the oxygen octahedra surrounding the Mn 
cations^ as well as a collective tilting of these octahe- 
dra, the so called "GdFe0 3 -type" distortion.^ The mag- 
netic moments of the Mn cations in LaMnC>3 order at 
Xn = 140 K in a so-called "A-type" antiferromagnetic 
structure P with parallel alignment of all moments within 
a certain (001) plane and antiparallel alignment of the 
moments between adjacent planes. 

In the ideal cubic perovskite structure (see Fig. |TJd) 
the 3d states of the Mn cations are split by the crystal 
field into the lower-lying three-fold degenerate t^g states 
and the higher-lying two-fold degenerate e g states. In 
LaMnC"3 the majority-spin i 2g states are fully occupied, 
whereas the two majority-spin e„ levels are filled with 




FIG. I: a) Experimentally observed Pnma structure of 
LaMnOs according to Ref. 7, b) Ideal cubic perovskite struc- 
ture. The oxygen anions form a network of corner-shared 
octahedra. The Mn cations (not shown) are situated in the 
centers of the oxygen octahedra and the La cations occupy 
the space between the octahedra. 



only one electron, according to the formal high-spin d 
electron configuration of the Mn 3+ cation. 

The theoretical modeling of manganite systems is usu- 
ally based on the assumption that the important low en- 
ergy dynamics of these systems can be described within 
an effective two band tight-binding (TB) model for the 
Mn e g electrons. In this model, electrons can hop be- 
tween the e g orbitals at neighboring Mn sites, and the 
corresponding kinetic energy term in the Hamiltonian is 
eventually complemented by local terms describing the 
coupling to the tig "core spin" , the coupling to the JT 
distortion of the surrounding oxygen octahedron, and the 
electron-electron interaction. These models can account 
for many of the properties observed in manganite sys- 
tems (see Ref. |5] for a recent review). Reasonable values 
for the model parameters, which describe the strengths 
of the various competing interactions, can be inferred 
from experiments. Nevertheless, it is very desirable to in- 
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dependently calculate values for these parameters based 
on fundamental "first principles" theory. Evaluating the 
models using such independently determined parameters 
provides a stringent test for the accuracy of the model 
description. 

Electronic structure calculations based on density 
functional theory (DFT^pEH provide a way to study 
the ground state electronic structure of a specific mate- 
rial without having to resort to model assumptions, and 
therefore provide a good starting point for the derivation 
of more simplified models and their parameters (see e.g. 
Refs . [TBI and 1X3)1 . The electronic structure of LaMnOs has 
been studied previously within the local spin density ap- 
proxim ation (LSD A) to DFT and by using the LSDA+C/ 
method? 15 * 16 * 17 * 18 -! It was shown that many properties 
such as the correct magnetic ground state and even some 
spectral properties are well described by these methods, 
provided the correct experimental crystal structure is 
used in the calculation. 

Although the model treatment of manganite systems 
usually employs a pure e g electron description, it is gen- 
erally understood that the electron hopping between the 
e g states on neighboring Mn sites is truly an effective 
hopping which is mediated by the intermediate oxygen 
anions via d-p or d-s hopping. The resulting bands with 
predominant e g character can be described by an effective 
two-band model if the Mn e g states are energetically sep- 
arated from the oxygen p and s states. In this case, the 
effective nearest neighbor hopping amplitude t between 
the Mn e g states is (to leading order) given by: 



t oc 



Ed — E„ 



(i) 



Here, t p d is the hopping amplitude between the Mn e g 
and the oxygen p states, Ed and E p are the energies of 
the corresponding ionic levels, and for simplicity we have 
neglected hopping via the oxygen s states. 

The JT distortion changes the Mn-0 bond lengths 
while the octahedral tilts change the bond angles; thus 
both distortions affect the overlap integrals which de- 
termine the hopping amplitude t p d- It is therefore not 
clear a priori that a simple effective TB model with 
fixed (distortion-independent) hopping amplitudes t can 
be used to study the effects of lattice distortions in man- 
ganite systems. 

Here, we use the Kohn-Sham band-structure calculated 
within the LSDA as a reference for the non-interacting 
TB model, and we analyze how well the relevant part 
of the energy dispersion of LaMn03 can be fitted within 
an effective two-band TB model for the e g electrons. In 
particular, we analyze the effects of the two dominant 
structural distortions in LaMnOs, the JT distortion and 
the GdFeC>3-type rotations, and we address the question 
of whether magnetic and orbital (JT) order affects the 
effective hopping amplitudes. 

The result of our analysis is that the effective two- 
band model gives a good fit of the e 9 -projected Kohn- 
Sham band-structure, provided that hopping between 



both nearest and next-nearest neighbors is taken into ac- 
count. We show that the same hopping amplitudes can 
be used for the ferromagnetic, the A- type antiferromag- 
netic, and the JT distorted case, so that the simple two- 
band TB model can be used to study the effects of JT dis- 
tortions. Furthermore we quantify the dependence of the 
hopping amplitudes on volume changes and on GdFeC>3- 
type rotations. The latter lead to significant reductions 
of the hopping amplitudes (~ 25 % for the experimental 
structure) relative to the ideal cubic structure with the 
same unit cell volume. The hopping amplitudes corre- 
sponding to the observed bond angles should therefore 
be used in theoretical modeling. 

Our results also provide a quantitative determination 
of the JT and Hund's rule couplings. The result for 
the Hund's coupling is consistent with previous work; 
the JT coupling is considerably smaller than previous 
estimatesPlMaij We find that both the JT and Hund's 
coupling are required to stabilize the insulating state 
within LSDA. 

Our conclusions rely in an essential way on the energy 
separation of the transition metal d-bands and the oxy- 
gen p-bands; methods such as LSDA+J7 which shift the 
energy of the transition-metal d-bands relative to the en- 
ergy of the oxygen p-bands can produce a band structure 
that is very poorly described by a simple two-band TB 
model. 

The remaining part of this paper is organized as fol- 
lows. Sec. [II] contains a brief summary of the methods 
and technical details of our work. We first describe the 
method we use for our LSDA calculations, then specify 
the TB Hamiltonian, and finally describe how we decom- 
pose the various structural distortions found experimen- 
tally in LaMnOs. Our results are discussed in Sec. HI 
and we end with a summary of our main conclusions and 
implications for future work. 



II. METHODS AND TECHNICAL DETAILS 
A. Computational method 

We calculate the LSDA Kohn-Sham band-structure 
for LaMnOs with both ferromagnetic and A-type anti- 
ferromagnetic order in various structural modifications 
using the projector augmented-wave (PAW) method im- 
plement edinthe "Vienna Ab-initio Simulation Package" 
(VASP)P^MIWe treat the La 5s, La bp, and Mn 3p 
pseudo-core states as valence states, and we use a plane- 
wave energy cutoff of 400 eV in all our calculations. We 
employ T-centered 6x6x6 and 4x4x3 fc-point grids for 
the calculations corresponding to the simple and quadru- 
pled perovskite unit cells, respectively, and correspond- 
ing grids for the structures in which the unit cells are 
doubled along the z direction or within the x-y plane. 
These values result in a good convergence of the calcu- 
lated band-structures. 

In order to extract the bands resulting from the Mn 
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e g states we use the "fatbands" technique, i.e. we assign 
a weight to each calculated eigenvalue, which is propor- 
tional to the amount of Mn e g character contained in 
the corresponding Bloch function, and we identify the 
e 9 -derived bands as those containing a non-negligible e g 
character. 

It has been shown in Ref. [TS] that the LSDA gives a 
good overall account of the electronic properties of man- 
ganite systems, even though the tendency to open up an 
energy gap between occupied and unoccupied states is 
underestimated within the LSDA. This is a well-known 
feature of the LSDA, which results from the inability of 
the LSDA to correctly account for the strong Coulomb 
correlations between the rather localized d states in tran- 
sition metal oxides. Such local Coulomb interactions are 
usually incorporated in the model Hamiltonian via a sep- 
arate interaction term. In the following we do not include 
such an interaction term in our model analysis, and thus 
the corresponding deficiencies of the LSDA do not af- 
fect our results (assuming that the separate treatment 
of local correlations is justified), except for the question 
related to the energy separation between the Mn d and 



the oxygen p states, which is discussed in Sec. IV 



B. Model Hamiltonian 

In Sec. IIIII we relate the calculated LSDA band- 
structure to the following TB model, which contains the 
terms that are typically used for the theoretical modeling 
of manganite systems (see e.g. Ref. [5]): 



H = H\ 



kin 



Hund 



H. 



JT 



with 
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(2) 

(3) 
(4) 

(5) 



Here, dg a a is the annihilation operator for an e g electron 

at site R in orbital a with spin a, and the correspond- 
ing boldface symbol indicates a pseudo-spinor in orbital 
space dig = (dg 1 a ,dg 2 a ) T . The orbital indexes 1 and 
2 correspond to |3z 2 — r 2 ) and \x 2 — y 2 ) orbitals, respec- 
tively. t^ + g = X^ =0 t^ + f l are the hopping ampli- 
tudes between site R and R+5 and r J are the usual Pauli 
matrices supplemented by the 2x2 unit matrix. Sj^ is the 

normalized core spin of the ti g electrons (|<5s| = 1), and 
Q x ^ are the amplitudes of the two JT modes at site R 





FIG. 2: Schematic depiction of the JT modes Q x (left side) 
and Q z (right side). The displacements Xg, Yg, and used 
in Eqs. |6f and {?} are indicated. Since we are considering 
only inversion symmetric distortions of the oxygen octahedra 
it is enough to specify the displacements of the oxygen anions 
adjacent to the central Mn cation in the positive x, y, and 
z directions. Note that |X^| = |Yg| in the case of Q x and 
for Q\ 



\Yn\ 



Zr\ 



that couple to the e g electrons: 
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(6) 
(7) 



Here, Xg, Yg, and are the displacements along x, y, 
and z of the oxygen anions that are situated adjacent to 
the Mn site at R in x, y, and z direction, respectively, 
and only inversion symmetric distortions of the oxygen 
octahedra are taken into account (see Fig. [2|. J and A 
are coupling constants for the local interaction terms. 

The first term in Eq. ^ describes the hopping between 
neighboring Mn sites. We will consider hopping between 
both nearest and next nearest neighbors. Symmetry dic- 
tates that the hopping matrices for nearest neighbor hop- 
ping are: 
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and for next nearest neighbor hopping: 
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Here, a is the lattice constant of the underlying cubic 
perovskite lattice. 
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FIG. 3: Staggered Q^-type JT distortion within the x-y plane 
found experimentally in LaMn03. The long and short axes 
of the octahedra alternate along the x and y directions, as 
shown in b), compared to the ideal structure shown in a). 

The second and third terms in Eq. Q describe the 
Hund's rule coupling to the ti g core spin and the coupling 
to the JT distortion of the oxygen octahedra surround- 
ing site R, respectively. The normalized ti g core spin is 
treated as a classical vector which is fixed externally in 
our model calculations. The values of Q X J Z are given by 
the positions of the oxygen anions used in our LSDA cal- 
culations and are also treated as external parameters in 
the TB model. 



C. Structural decomposition 

As described in the introduction, LaMnC>3 exhibits a 
strongly distorted perovskite structure with space group 
PnmaW The deviation of this structure from the perfect 
cubic perovskite structure (with space group Pm3m) can 
be decomposed into the following three contributions: 

(i) A staggered (checkerboard-like) Q^-type JT distor- 
tion of the oxygen octahedra within the x-y plane, 
with the long and short axes of neighboring octa- 
hedra alternating between the x and y directions 
(see Fig. |3| . This JT distortion leads to a doubling 
of the unit cell compared to the ideal cubic per- 
ovskite structure, with new "in- plane" lattice vec- 
tors a = ao(x — y) and b = ao(x+y), where ao is the 
lattice constant of the original (undistorted) per- 
ovskite structure. Identical x-y planes are stacked 
on top of each other along the z direction. The 
resulting symmetry is tetragonal. 

(ii) "GdFeC>3-type" rotations (tilting) of the oxygen oc- 
tahedra, leading to an additional doubling of the 
unit cell along the z direction, with the new lattice 
vector c = 2a^z, and a reduction to orthorhombic 
Pnma symmetry. 

(iii) Displacements of the La cations, and a deformation 
(strain) of the parallelepiped formed by the lattice 



TABLE I: Decomposition of the experimentally observed 
distortion of the oxygen network into JT component and 
GdFeOs-type (GFO) rotations. The upper four lines contain 
the Wyckoff positions of the two inequivalent oxygen sites Ol 
(4c) and 02 (8d) in the experimental structure (Ref. [7)), in 
the ideal cubic perovskite structure (note that these values 
correspond to a quadrupled unit cell), and our decomposition 
in pure JT and pure GdFe03-type components. Ar, ArjT, 
and Ar*GFo represent the corresponding full experimental dis- 
tortion, and its decomposition into pure JT and GFO-type 
distortion, respectively, x, y, and z are the coordinates with 
respect to the orthorhombic lattice vectors. 
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0.25 
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0.0385 


Ar 


-0.0733 


-0.0107 


0.0 


-0.0243 


0.0514 


0.0385 


AfjT 


0.0 


0.0 


0.0 


0.0136 


0.0136 


0.0 


AFgfo 


-0.0733 


-0.0107 


0.0 


-0.0379 


0.0379 


0.0385 



vectors a, b, and c, consistent with the orthorhom- 
bic crystal class. 

We expect that the internal distortions of the oxygen 
network, i.e. components (i) and (ii) described above, 
have the largest effect on the e g bands of LaMnC>3 via the 
ligand-field splitting, whereas the influence of the lattice 
strain and of the La displacements, i.e. component (iii), 
can be neglected. We test the validity of this hypothesis 
in Sec. ESI 

In the following, except for the test calculations pre- 
sented in Sec. |III A[ we therefore consider only compo- 
nents (i) and (ii). This means that we use the "pseudo- 
cubic" lattice vectors a, 6, and c defined above, discard- 
ing any orthorhombic strain, and we place the La cations 
on their ideal positions, corresponding to the perfect cu- 
bic perovskite structure. For the internal distortion of 
the oxygen network, we use the experimental data ob- 
tained in Ref. which we decompose into the pure JT 
distortion and the GdFeC>3-type distortion, as described 
below. Furthermore, we use a "cubic" lattice constant ao 
= 3.9345 A, which results in the same volume per formula 
unit as in the experimentally observed structured 

Table |l] lists the Wyckoff positions for the two inequiv- 
alent oxygen sites Ol (4c) and 02 (8d) in the experimen- 
tally determined Pnma structure, and how we decom- 
pose the corresponding structural distortion in the pure 
JT component (i) and the GdFe03-type distortion (ii). 
The decomposition is such that Ar = Ar*jx + Atqpq 
and ArjT is orthogonal to Aropo, where Ar, AfjT, and 
ArcFO are the full experimental distortion and its de- 
composition into pure JT and GdFe03-type distortion, 
respectively. Since the Wyckoff coordinates x, y, and 
z can be directly interpreted as the coordinates relative 
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to the lattice vectors a, b, and c, it follows that in the 
purely JT distorted structure each oxygen anion 02 is 
displaced by |Ar JT (02)| = |Az JT (02) a + Ay JT (02) b + 
Az JT (02) c I = 0.1070 A. According to Eq. (j6) this corre- 
sponds to a JT amplitude of Qg = 0.1513A = 0.0385 do- 



III. DISCUSSION OF RESULTS 

A. Test of structural decomposition 

In the previous section we stated that only components 
(i) and (ii), i.e. the internal distortion of the oxygen net- 
work, are important for the e g bands in LaMnOa, and 
that the lattice strain as well as the displacements of the 
La cations are negligible. In order to test this hypothe- 
sis, we now compare the LSDA band-structure calculated 
for the full experimental structure of Ref. [7] with the one 
calculated for the slightly simplified structure described 
above, where the lattice strain and the La displacements 
are set to zero, while the internal coordinates of the oxy- 
gen anions are the same as observed experimentally. 

The corresponding LSDA band-structures in the en- 
ergy range of the Mn e g bands calculated for A-type an- 
tiferromagnetic ordering are shown in Fig.[4]along certain 
high symmetry directions of the orthorhombic Brillouin 
zone. The Mn e g bands are visualized by dots along 
the bands, with the radius of the dots proportional to 
the amount of Mn e g character contained in the corre- 
sponding Bloch-function. It is clearly seen that the band- 
structures obtained for the fully experimental structure 
and for the simplified structure with only the oxygen dis- 
tortions included are nearly indistinguishable, with only 
small deformations of the energy bands resulting from 
the orthorhombic strain and the La displacements. This 
validates our initial hypothesis, and in the following we 
therefore analyze only the effect of the internal structural 
distortion of the oxygen octahedra on the dispersion of 
the e g bands. 

We point out that by setting the lattice strain to zero 
we also neglect any homogeneous Q z -type JT distortion. 
The good agreement between the two band-structures 
shown in Fig. [4] thus also indicates that there is no no- 
ticeable effect of Q z on the electronic band-structure of 
LaMnC>3 in its experimental crystal structure. 

The simplified structure that gives rise to the LSDA 
band-structure shown in the bottom part of Fig. [4] re- 
sults from the superposition of distortions (i) and (ii) 
(described above) of the oxygen network. In the following 
we will first establish the e g band-structure of LaMnC>3 in 
the ideal cubic perovskite structure and then separately 
introduce either the JT distortions, component (i), or the 
GdFeC>3-type rotations, component (ii), and analyze the 
corresponding effects on the e g bands. 
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FIG. 4: LSDA energy bands (thin lines) along high symme- 
try directions of the orthorhombic Brillouin zone calculated 
for a) the exact experimental structure, and b) for the simpli- 
fied structure with no orthorhombic strain and the La cations 
on their ideal positions. Both calculations are done for A- 
type antiferromagnetic ordering. The thick dots indicate the 
amount of Mn e 9 character in the corresponding Bloch func- 
tions. The dashed horizontal line at ~ 7.75 eV indicates the 
Fermi energy. In b) the high symmetry fc-points are given 
in cartesian coordinates and in units of n/ao, in a) the cor- 
responding fc-points are labeled identically but correspond 
to the slightly strained reciprocal lattice of the experimen- 
tal structure. The two cases are nearly indistinguishable. 



B. Cubic structure 

Fig. [5] shows the calculated LSDA energy dispersion 
around the Fermi-level, calculated for ferromagnetic cu- 
bic LaMnC>3 in the undistorted cubic perovskite struc- 
ture with ao = 3.9345 A. Only the majority spin bands 
are shown. Within the TB model defined in Eq. Q, 
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FIG. 5: Majority spin bands for cubic ferromagnetic LaMnC>3. 
The LSDA band-structure is represented by thin lines, with 
the dots indicating the amount of Mn e 9 character in the 
corresponding Bloch function. The thick broken line is the TB 
fit with only nearest neighbor hopping, whereas the thick solid 
line represents the TB fit including both nearest and next 
nearest neighbor hopping. The thin dashed horizontal line 
at ~ 8.2 eV indicates the Fermi energy. The high symmetry 
k-points are given in cartesian coordinates and in units of 
7r/ao. 



the Hund's coupling simply splits the spin majority and 
minority bands rigidly by AE = 2 J for a ferromagnetic 
arrangement of the ti g core spins, and we therefore dis- 
cuss only the majority spin bands in the following. The 
Mn e g bands are again visualized by the dots along the 
bands in Fig. [5] It is evident that even though the Mn e g 
bands are intersected by other bands, the e g dispersion 
can be nicely traced along the dots. The e g bands are 
about half-filled, as expected from the formal electron 
configuration t\ g e g of the Mn 3+ ion. The bands at ~ 
7eV and lower, just touching the lower range of the e g 
bands, are the filled majority Mn t2 g bands. The weakly 
dispersive bands slightly below 10 eV that intersect the 
Mn e g bands correspond to the La 4/ states, and the 
strongly dispersive unoccupied bands above the Mn e g 
manifold have predominantly La d character. 

The thick lines in Fig. [5] correspond to fits of the 
nearest and next nearest neighbor TB models for the 
e g bands. The nearest neighbor hopping parameter 
t = 0.655 eV is determined from the full e g bandwidth 
W = 3.928 eV= 6t. The next nearest neighbor hopping 
parameter t' is obtained in the following way: In the next 
nearest neighbor model the width of the energy disper- 
sion of the upper e g band between /c-points T = (0, 0, 0) 

(2) 

and X = (1, 0, 0) is equal to AE^ — At — 16i', whereas 
the dispersion width of the lower band between X and 
M = (1,1,0) is equal to AE^ M = 2t- 16t': 25 The corre- 
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FIG. 6: Band-structure for cubic LaMnOa with A-type anti- 
ferromagnetic order. The LSDA band-structure is represented 
by thin lines, with the dots indicating the amount of Mn e 9 
character in the corresponding Bloch function. The thick line 
represents the TB fit including both nearest and next nearest 
neighbor hopping. The thin dashed horizontal line at ~ 8.3 eV 
indicates the Fermi energy. The high symmetry k-points are 
given in cartesian coordinates and in units of 7r/ao. 



sponding energy differences obtained from the LSDA cal- 
culation are AE^ = 0.851 eV and AE^ M = 1.965 eV, 
leading to t' = 0.041 eV or t' = 0.029 eV, respectively 
(and using the previously obtained £=0.655 eV). An av- 
erage value of t! — 0.035 eV is used for the TB fit in Fig. [5] 
and in the remaining part of this paper. 

It becomes clear from Fig. [5] that the simple nearest 
neighbor TB model cannot reproduce the LSDA disper- 
sion very well, whereas the next nearest neighbor TB 
model leads to a very good description of the energy dis- 
persion for all fc-points. 

We point out that a nonmagnetic LDA calculation re- 
sults in a low-spin electron configuration and the loss of 
the ti g core spin. Indeed, if we perform a nonmagnetic 
LDA calculation, the e g bands are empty and higher in 
energy compared to the oxygen p levels, which in accor- 
dance with Eq. results in a reduced bandwidth of 
3.546 eV, corresponding to a nearest neighbor hopping 
amplitude of t = 0.591 eV. A nonmagnetic LDA calcu- 
lation is thus not necessarily a good representation of 
the electronic structure of the paramagnetic phase, and 
we therefore use the ferromagnetic state as the starting 
point for the model analysis. In general, this shows that 
an LDA+DMFT treatment of LaMnC>3 based on a non- 
magnetic LDA calculation, such as the one presented in 
Ref. |571 leads to a slight underestimation of the electron 
hopping. 

Next we investigate the influence of A-type antiferro- 
magnetic order. Fig. [6] shows the calculated LSDA band- 
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structure for this case. Note that the underlying crystal 
structure is still perfect cubic perovskite. Again, the e g 
character of the bands is visualized by the dots. The 
thick lines corresponds to the fit within our antiferro- 
magnetic next nearest neighbor TB model with the hop- 
ping parameters obtained from the ferromagnetic case. 
Due to the doubling of the magnetic unit cell the num- 
ber of bands is also doubled. A Hund's-rule parameter 
,/ = 1.340 eV is obtained from the energy splitting at 
the r point between the two bands at 6.7 eV and 9.4 eV, 
which show no dispersion along TA (A = (0,0,1/2)). 
This splitting is exactly equal to 2J in the TB model. 
The value J — 1.340 eV is within 2.4% of the value 
J = 1.309 eV obtained from the ferromagnetic band- 
structure as the splitting between majority and minority 
spin states at the T point (not shown). 

Fig. IBIshows that the two lowest e g bands are described 
very well by the antiferromagnetic TB model. The upper 
two bands show some deviations from the model, espe- 
cially in the high energy region. This is an inevitable 
result of the description within a pure d band model. As 
described in Sec. |TJ the "true" hopping is mediated by 
the oxygen p orbitals and therefore the e g dispersion de- 
pends on the energetic distance from the oxygen p levels 
(see Eq. 0). This leads to a slight overestimation of the 
energy dispersion for the high energy states in the pure 
d model. The same effect can also be observed in the 
ferromagnetic case: Due to their higher energy relative 
to the oxygen p states, the bandwidth of the e g minority 
spin bands is smaller than for the corresponding major- 
ity spin bands. The nearest neighbor hopping parameter 
corresponding to the minority spin bands in the ferro- 
magnetic case is t = 0.548 eV. In the following we use 
the value t = 0.655 eV, corresponding to the majority 
spin bands in the ferromagnetic configuration, since this 
value is representative for the e g bands close to the Fermi 
level which determine the important low energy behavior 
in manganite systems. 

Fig. [7] shows the dependence of both nearest and next 
nearest hopping parameters on the lattice constant a. 
The nearest neighbor hopping t decreases with increas- 
ing Mn-0 bond length, whereas the next nearest neigh- 
bor hopping t' shows a slight increase. This somewhat 
unexpected behavior of t' results from the fact that the 
energy difference between oxygen p and Mn d states de- 
creases with increasing volume, and therefore counteracts 
the effect of the reduced overlap integrals for larger lat- 
tice constants. 



C. Purely JT distorted structure 



We now address the effect of the JT distortion of the 
oxygen octahedra on the e g bands in LaMnC>3. As de- 




scribed in Sec. II C we consider only the internal dis- 
tortions of the oxygen network. This means that we 
displace the oxygen anions relative to the cubic struc- 
ture according to the decomposition of the Wyckoff posi- 



FIG. 7: Variation of hopping amplitudes with lattice constant. 
Open circles correspond to the nearest neighbor hopping t 
and filled squares correspond to the next nearest neighbor 
hopping t' . a indicates the lattice constant derived from the 
experimental unit cell volume. 



tions described in Table [I] (line "JT"). This results in the 
staggered arrangement of JT distorted oxygen octahedra 
shown in Fig. J3Jd, with the same JT amplitude Qq as in 
the experimental structure. 

Fig. [8] shows the calculated LSDA majority spin bands 
for LaMnC>3 in the purely JT distorted structure with fer- 
romagnetic spin ordering. Due to the unit cell doubling 
within the x-y plane compared to the cubic structure, 
the two bands corresponding to the fc-points between 
X — (1, 0, 0) and M = (1, 1, 0) in the ferromagnetic cubic 
case (see Fig. [5]) are "back-folded" between the /c-points 
U = (1, 0, 0) and T = (0, 0, 0) in the tetragonal Brillouin 
zone. In addition, the two bands between M — (1,1,0) 
and R = (1,1,1) in Fig. plnow correspond to the two up- 
per bands between r = (0, 0, 0) and Z = (0, 0, 1) in FigM 
and the two bands between U = (1, 0, 0) and R = (1,0,1) 
in Fig. [8] are now twofold degenerate. It can be seen that 
the level splitting between the two original e g states at 
the r point (at ~6.7eV in Fig. [8]), which are degenerate 
in the cubic case, is very small (~ 0.08 eV), and that the 
main effect of the JT distortion is to remove the band 
crossing between T and U — (1,0,0), resulting from the 
simple "back-folding" of the cubic band-structure due to 
the unit cell doubling. 

To obtain the value of the JT coupling constant A 
within our TB model, we first determine the fc-point of 
the band crossing between T and U for the case of zero 
JT distortion. We then determine XQq in the model by 
fitting the splitting at this fc-point to the corresponding 
splitting obtained from the LSDA calculation. In this 
way we obtain a value of XQq — 0.248 eV, correspo ndin g 
to A = 1.639 eV/A, since Qg = 0.1513 A (see Sec. |ll C|. 
Alternatively we can also fit the small splitting of the two 



lowest e g bands at the T point by numerically adjusting 
the JT coupling. In this way we find XQq = 0.289 eV, 



(0,0,1) (0,0,0) (1,0,0) (1,0,1) 

FIG. 8: Majority-spin band-structure for ferromagnetic 
LaMn03 in the purely JT distorted structure (see Fig. [3Jd) , 
where the unit cell is doubled in the x-y-p\a,ne. The high 
symmetry k-points are given with respect to the cartesian co- 
ordinate system defined by the reciprocal lattice vectors of the 
undistorted cubic structure. The units are n/ao- The LSDA 
band-structure is represented by thin lines with the dots in- 
dicating the amount of Mn e g character in the corresponding 
Bloch function. The thick line represents the TB fit and the 
thin dashed horizontal line at ~ 8.1 eV indicates the Fermi 
energy. 

corresponding to A — 1.910 eV/ A, which is within 17% 
of the value obtained above. This shows that the ex- 
tracted coupling strength does not depend critically on 
the fitting procedure. 

The energy dispersion calculated within the TB model 
using the hopping amplitudes t — 0.655 eV and t' = 
0.035 eV obtained for the cubic structure and the JT 
coupling constant A = 1.639 eV/ A obtained in the way 
described above is compared to the full LSDA band- 
structure in Fig. [8j It can be seen that the dispersion of 
the e g -projected bands is well described within the TB 
model. Some deviations occur close to the lifted band 
crossing between V and U — (1,0,0). These deviations 
are most likely caused by the asymmetry of the Mn-0 
bonds, which is neglected in the effective e g TB model. 
The good overall quality of the TB fit shows that the TB 
hopping amplitudes are not affected by the presence of 
the JT distortion. This indicates that the model descrip- 
tion with the assumed local coupling of the JT distortion 
to the e g levels is justified. 

To test whether the linear form of the coupling term 
within our model is consistent with the density functional 
theory calculation, and to further test our procedure for 
obtaining XQ X from the LSDA results, we perform ad- 
ditional calculations with different amplitudes of the JT 
distortion. The results are presented in Fig. [9j where the 
JT energy XQ X is determined by fitting the band-splitting 



FIG. 9: Magnitude of XQ X obtained by fitting the TB model 
to LSDA calculations with different amplitudes of the JT dis- 
tortion. Q x indicates the magnitude of the JT distortion 
found in the experimental structure. Filled circles are the 
results of the actual calculations, whereas the straight line 
corresponds to XQ X = 0.247 eV Q x /Q x . 



between V and U — (1,0,0). We note that XQ X on the 
ordinate in Fig. [9] should be regarded as the model pa- 
rameter that we obtain by our fitting procedure, whereas 
the abscissa Q x /Q x characterizes the input structure for 
our LSDA calculation (relative to the experimentally ob- 
served JT distortion Q x ). It is evident that the depen- 
dence of XQ X on the input distortion is nearly perfectly 
linear, which indicates the good quality of our fit and the 
adequacy of the linear coupling term within the model. 

It becomes clear from Fig. [8] that the JT distortion 
in the experimentally observed structure of LaMnC>3 has 
only a weak effect on the dispersion along k z . This is 
complementary to the effect of the A-type magnetic or- 
der, which strongly suppresses the electron hopping along 
this direction. The insulating band-structure obtained in 
LSDA for the fully distorted structure with A-type an- 
tiferromagnetic order is therefore a combined effect of 
both the staggered JT distortion within the x-y plane 
and the A-type antiferromagnetism. To achieve an insu- 
lating state within LSDA solely due to the JT distortion 
(i.e. for the ferromagnetic case), would require an un- 
realistically large JT amplitude. Within our TB model, 
a value of XQ X > 1.1 eV, i.e. more than four times the 
JT distortion of the experimental structure, is required 
to open up an energy gap. This is due to the the large 
value of the hopping t and the fact that for staggered JT 
order Hjt does not commute with H^ n . The fact that 
the JT distortion alone is not enough to stabilize an in- 
sulating state in LaMnC>3 has also been pointed out in 
Refs. I19I21I28I 

One possibility for explaining the insulating charac- 
ter of LaMnC>3, as noted by previous authors™ is that 
electron-electron interactions beyond LSDA increase the 
effective JT splitting, thereby stabilizing the insulating 



9 




(0,0,1) (0,0,0) (1,0,0) (1,0,1) 



FIG. 10: Majority-spin band-structure for ferromagnetic 
LaMn03 in the purely JT distorted structure calculated us- 
ing the LSDA+f/ method with U c « = 7.12 eV (U = 8eV and 
J = 0.88 eV). The high symmetry fc-points are given in 
cartesian coordinates and in units of n/ao. The Kohn-Sham 
band-structure is represented by thin lines with the dots in- 
dicating the amount of Mn e g character in the corresponding 
Bloch function. The thin dashed horizontal line at ~7.2eV 
indicates the Fermi energy. The overlapping bands in the 
energy range between eV and 6 eV have mixed O 2p-Mn 3d 
character. 
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(0,0,1/2) (0,0,0) (1,0,0) (1,0,1/2) 

FIG. 11: LSDA energy bands (majority spin) for ferromag- 
netic LaMnOa with the experimentally observed GdFeOs- 
type tilting of the oxygen octahedra, but without JT distor- 
tion. The LSDA band-structure is represented by thin lines 
with the dots indicating the amount of Mn e g character in the 
corresponding Bloch function. The thick line represents the 
TB fit with reduced hopping amplitudes and the thin dashed 
horizontal line at ~ 7.8 eV indicates the Fermi energy. The 
high symmetry k-points are given in cartesian coordinates and 
in units of n/ao- AEi and AE2 mark the energy differences 
plotted in Fig. [12] 



state. To address this we have performed additional 
LSDA+f/ calculations (which will be discussed in detail 
elsewhere). In the LSDA+[/ method, the interactions 
between the d states of the transition metal cations are 
treated explicitly on the level of a mean-field Hubbard 
modelP^ Thereby, a parameter U represents the strength 
of the (screened) on-site Coulomb repulsion between the 
d electrons and a parameter J H represents the Hund's 
coupling. In our LSDA+J7 calculations we use a slightly 
simplified approach where only [/ e fr = U — J H enters.^ 
We expect that the on-site Coulomb repulsion en- 
hances the effect of the JT distortion and therefore drives 
the system towards an insulating state as the value of U 
is increased. However, the calculated LSDA+t/ band- 
structure for the purely JT distorted case with Q x = Qq 
and ferromagnetic spin order stays metallic even for 
a rather large Hubbard parameter of U = 8eV (see 
Fig. 10 1. The reason for this is the following: The use of 



the LSDA+C/ method leads to a strong downward energy 
shift of the occupied part of the e g bands and also to a 
moderate overall downshift of the d states as a whole. 
As a result of the latter, the higher-lying, mostly unoc- 
cupied e g bands in LaMn03 move somewhat closer to 
the O 2p bands, which are located in the energy range 
between OeV and 6eV and are not shifted within the 
LSDA+J7 method. Furthermore, the strong downward 
shift of the occupied d bands results in an energetic over- 



lap and therefore strong hybridization between the oc- 
cupied e g states and the O 2p bands. In this case the 
simple two-band TB model is not applicable any more, 
and some mixed p-d bands, which extend above the orig- 
inal top of the O 2p states at 6 eV, reach above the Fermi 
level. Since the LSDA+t/ method corrects only for the 
Coulomb interactions between the Mn d states, it is not 
obvious whether this shift of the d states relative to the 
O 2p states is a real physical effect or rather an artifact 



of the LSDA+[/ method (see Sec. IV for a further discus 



sion of this point). An exact experimental determination 
of the energy separation between the Mn d and Op states 
in LaMnOa would provide further insight on this. 



D. GdFeC>3-type rotations 

Finally, we address the effect of the GdFe03-type oxy- 
gen octahedra rotations. These distortions change the 
Mn-O-Mn bond angles and are therefore expected to al- 
ter the magnitude of the hopping amplitudes. In addi- 
tion, due to the resulting symmetry lowering, this dis- 
tortion will enable hopping between orbitals that was ei- 
ther symmetry-forbidden or negligibly small in the undis- 
torted state. 

Fig. [TT] shows the calculated LSDA energy bands for 
the structure where only the GdFe03-type distortion, 



10 




170 160 150 
Mn-O-Mn bond angle [degrees] 



FIG. 12: Reduction of some characteristic energy differences 
at the fc-point U = (1,0,0) (see Fig. Ill for different ampli- 
tudes of the GdFe03-type distortion, which are indicated by 



the corresponding Mn-O-Mn bond angle. 



component (ii) in Sec. II C is included, whereas the JT 
distortion, component (i), is set to zero. The Wyckoff 
positions of the oxygen anions for this configuration are 
listed in Table [I] (line "GFO"). In this structure the Mn- 
O-Mn bond angles are reduced from the ideal 180° to 
about 155°. The thick line in Fig. 11 corresponds to a 



fit within the next-nearest neighbor TB model with both 
nearest and next-nearest hopping amplitudes scaled by a 
factor of 0.7 compared to the perfectly cubic case. The 
high quality of the fit is striking, even though the LSDA 
band structure shows some additional dispersion at the 
bottom of the e g bands which is not accounted for in 
the TB model. It appears that, to a good accuracy, the 
oxygen tilts can be incorporated in the model simply by 
reducing the hopping amplitudes in an appropriate way 
without having to include additional hopping parameters 
due to the lower symmetry. 

To further quantify the reduction of the hopping am- 
plitudes as a result of the GdFeC>3-type distortion, we 
perform calculations for different degrees of distortion, 
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by scaling Ar*QFO (see Table [I]) accordingly. Fig 
shows the resulting reduction of some characteristic en 
ergy differences at the fc-point U = (1, 0, 0) as a function 
of the Mn-O-Mn bond angle. If the GdFe03-type rota- 
tions would lead to a simple scaling of the undistorted 



band-structure, all the lines in Fig. 12 would fall on top 
of each other. It can be seen that this is in fact a good ap- 
proximation for Mn-O-Mn bond angles down to ^155°, 
which corresponds to the experimentally observed struc- 
ture of LaMnOs. For bond angles smaller than 155° the 
band-structure starts deviating more significantly from 
the cubic case due to new hopping paths that become 
allowed in the distorted structure. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have shown, by comparing LSDA 
band-structure calculations to simple TB models, that 
the relevant electronic states in LaMnOs are well de- 
scribed by a model of e g orbit als with nearest and next- 
nearest neighbor hoppings. We have quantified the effect 
of changes in bond length (Fig. [7]) and of the octahe- 
dral rotations (Fig. 12 1 on the hopping parameters, and 
we find that for physically relevant values (bond angles 
> 155°) the GdFe03-type rotations significantly change 
the value of the hopping parameters but do not invalidate 
the TB description. Of particular importance is our find- 
ing that both the JT lattice distortions and the Hund's 
rule coupling are quantitatively accounted for by adding 
on-site interactions to the TB model, without changing 
the hopping parameters. In summary, these results jus- 
tify a TB plus interaction description of manganite sys- 
tems, and suggests more generally that such a description 
is useful for transition metal oxides. 

The parameters for nearest and next nearest neigh- 
bor hopping (defined in Eqs. (|3| and (|8|-(13l) which 
we obtain for the ideal cubic structure with lattice con- 
stant a = 3.935 A are t = 0.655 eV and t' = 0.035 eV. 
The Hund's rule and JT coupling constants (defined 
in Eqs. @-(|5])) which we obtain from our analysis are 
J = 1.34eVandA = 1.64eV/A. Our value of J = 1.34eV 
for the Hund's coupling is in excellent agreement with the 
exchange splitting 2 J = 2.7 eV derived from x-ray ab- 
sorption measure ments,^ and agrees well with previous 
LSDA calculations |I212Zl On the other hand, the strength 
of the JT coupling A obtained in this work is smaller than 
what has been concluded from x-ray absorption measure- 
ments in Ref. |3T] It is also smaller compared to various 
other values obtained previously from (mostly less elab- 
orate) fits of similar TB m odels to L SDA or LSDA+C/ 
band-structure calculations! 19 1 20 1 21 ^ 

Popovic and Satpathy used a fitting procedure very 
similar to the one presented here, and obtained a JT cou- 
pling strength of A = 2.8eV/A ( in our notation), com- 
pared to A = 1.64 eV/A obtained in the present work!^ 
We ascribe this difference to the use of the atomic sphere 
approximation (ASA) in the LSDA calculation of Ref. 20 
In the ASA the crystal volume is represented by an ar- 
rangement of overlapping "atomic spheres" This over- 
lap introduces an error, which furthermore depends on 
the amplitude of the JT distortion (since the JT distor- 
tion changes the overlap between the atomic spheres) and 
thus can have a pronounced effect on A. 

Ahn and Millis used a TB model very similar to the 
one used in this work, except that they didn't include the 
effect of next-nearest neighbor hopping.^ They obtained 
a value of A = 3.38 eV/A by simultaneously fitting 15 en- 
ergies at 4 different high symmetry /c-points to a previous 
LSDA calculation for the fully distorted antiferromagnet- 
ically ordered case. It is not obvious how sensitive such 
a simultaneous root mean square fit of all the model pa- 
rameters is to the exact value of A, but we expect that 
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the neglect of next nearest neighbor hopping will lead 
to a renormalization of the other parameters of the TB 
model in order to account for the missing dispersion due 
to the next-nearest neighbor hopping. 

A value of A = 2.85 eV/A was obtained by Yin et al.^ 
by calculating the dependence of several quantities on 
the amplitude of the JT distortion. This dependence 
was first obtained from LSDA+C/ calculations within a 
Wannier function representation, and then compared to 
the corresponding results calculated for a model Hamil- 
tonian including electron-electron interactions within the 
Hartree-Fock approximation. As in the case of Ref. [TH] 
discussed above, it is not clear how sensitive this simul- 
taneous fit of all parameters in the model Hamiltonian is 
to moderate changes in A. 

In contrast, the fitting procedure described in this 
work isolates the effect of each term in the Hamiltonian 
(Eq. Q) and thus allows to obtain each parameter inde- 
pendently from all others. In particular, it becomes clear 
from our calculated band-structure shown in Fig. [8] that 
the JT distortion does not lead to a rigid splitting of the 
e g bands, but that instead it has only subtle, albeit rather 
important effects on band crossings at certain fc-points 
in the Brillouin zone. The reason for this is that the JT 
Hamiltonian Hjt in Eq. ([5| for Q x =^ and staggered 
order does not commute with the kinetic energy term 
in Eq. (J3j> , and that due to the relatively small value 
XQq ss 0.25 eV the hopping energies are dominant, so 
that in general the effect of the JT distortion is only vis- 
ible as second order shifts in the energy. In other words, 
at a generic Ai-point the states picked out by the hopping 
term are not the eigenstates of the JT distortion. This 
suggests that the straightforward interpretation of peak 
splittings in the x-ray absorption spectra of Ref. |3T] as a 
direct consequence of the JT distortion is not necessarily 
justified. 

Finally, our analysis enables us to clearly identify the 
limitations of the effective two band e g TB description 
of manganite systems. Our TB analysis was success- 
ful because in LaMnC>3 within LSDA the e g bands are 
well-separated from the oxygen 2p bands, and neither 
the JT distortion nor the magnetic order change this en- 
ergy spacing, and thus the value of the effective hop- 
ping, significantly. The dependence of the effective hop- 



ping parameters on the energetic distance between the 
Mn e g and the O 2p states (see Eq. ([I])) is visible as 
a 15-20 % difference between the majority-spin and the 
minority-spin bandwidths and dispersion in the ferro- 
magnetic LSDA calculation for the cubic structure (see 
Sec. IIIBl, and also in the high- lying bands of the an- 
tiferromagnetic LSDA band-structure shown in Fig. [6] 
However, in cases where the e g and O 2p bands overlap in 
energy, such as for example in our LSDA+C/ calculation 
for the purely JT distorted structure shown in Fig. |10| 
the effective e g TB analysis fails, and the O 2p levels have 
to be taken into account explicitly. 

The energy shift of the occupied d states relative to 
the oxygen p states within the LSDA+C/ method is 
mainly caused by the so-called "double-counting correc- 
tion" , which attempts to correct for those contributions 
of the electron-electron interaction that are accounted for 
in both the LSDA and the local Hartree-Fock ("+£/") 
treatment. Since the double-counting correction is noto- 
riously ill-defined, this raises the question of whether such 
level shifts due to the electron-electron interaction and 
the resulting substantial renormalization of the effective 
hopping parameters are real effects, or whether this is an 
artifact of the LSDA+C/ scheme, which only accounts for 
the static (mean-field) electron-electron interaction be- 
tween the transition metal d states, while leaving the O 
2p states unchanged (we point out that the same prob- 
lem is also present within an LDA+DMFT treatment of 
electronic correlations). Optical evidence (see Ref. I54")l 
suggests that the O 2p bands in manganites are located 
about 4eV below the Fermi level, consistent with the 
LSDA result, but more detailed investigations of the en- 
ergy separation between the Mn e g and O 2p bands will 
be useful for future studies. 
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